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Abstract. We investigate the influence of a laser field, which is near-resonant to the atomic sodium 
3 2 Pi/2 — * 3 2 D3/2 transition, on the last bound levels of the A 1 state in Na2. In a molecular beam 
experiment level shifts up to m 100 MHz and light induced line broadenings were observed using an opti- 
cal double resonance excitation scheme. Moreover, the coupling laser can reduce the number of bound levels 
of the A state by one or more units, which effectively means that in the picture of a collision of a 3 2 Si/2 
and a 3 2 Pi/2-atom the scattering phase is altered by more than tv. The observed effects are interpreted as 
light induced couplings of the A 1 state, which correlates to the 3si/2+3pi/2 asymptote, to the 4 1 X'j~ 
and 2 1 n g states at the 3s!/2+3d3/2 asymptote. We performed multi-channel calculations, applying the 
mapped Fourier grid method, which reproduce our experimentally observed level shifts well. 
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1 Introduction 

Due to the fast progress in cooling and trapping of ul- 
tracold atoms and the achievement of Bose-Einstein con- 
densation 12,3 4,5,6,7], the interest in a detailed knowl- 
edge of ultracold collisional properties has increased [8]. 
They are often described by the s-wave scattering length 
a s , which is correlated to the scattering wave function 
at vanishing collisional energy. It has been shown [MTUl 
ITT] that this important physical quantity can be derived 
from spectroscopy of the vibrational levels and resonances 
around the atomic asymptote. But the scattering length is 
difficult to derive purely on theoretical grounds because it 
does not only depend on the long-range part of the molec- 
ular ground state interaction potentials but also on the 
accumulated phase from the inner part of the potential. 

Beside the knowledge of cold collision properties, the 
possibilities of their manipulation are central points of in- 
terest. It has been shown in Na-Na, Rb-Rb, and recently 
also in Cs-Cs ultracold collisions that the sign of the scat- 
tering length can be altered in the vicinity of a magnetic 
field induced Feshbach resonance by varying the magnetic 
field strength |12II13) . This offers a wide range of experi- 
ments on the dynamics in ultracold ensembles and BEC 
and a further understanding of two and three particle in- 
teraction. Unfortunately, magnetic tuning of two particle 
interaction is not applicable in typical magnetic traps be- 
cause the trapping field conditions are affected by the field 
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strength needed according to the desired tuning of the col- 
lision conditions. 

An alternative approach for the manipulation of cold 
collisions has been developed in ^1], where DC electric 
fields are proposed to control diatomic collisions at ul- 
tralow temperatures. 

A further promising proposal is the use of near reso- 
nant light fields to influence the scattering length a s |15l 
HrJinT) . This technique offers the possibility to manipulate 
the two-particle interaction by coherent light fields. First 
experiments with photoassociation |18| show that it is pos- 
sible to observe the optical analogue to the magnetic field 
induced Feshbach resonances if a laser with a frequency 
close to a photoassociation transition is focused into a 
magneto optical trap. The study of an optically induced 
Feshbach resonance is important for possible control of 
condensate dynamics. 

In this paper we follow another approach to determine 
the effect of near resonant laser light on cold atomic col- 
lisions. We use the high resolution of Doppler free spec- 
troscopy on a molecular beam to investigate the effect of 
a laser field that is near resonant to the atomic sodium 
3 2 Pi/2 ^3 2 D3/ 2 transition. Thus, this laser induces a cou- 
pling of the A 1 27+ state at the 3si/2+3pi/ 2 asymptote to 
states at the 3si/2+3d 3 / 2 asymptote (4 1 27+, 2 1 n g ) with 
their vibrational and continuum manifold. This results 
in energy shifts of asymptotic rovibrational energy levels 
with internuclear separations of about 100 A in all coupled 
states modeling atomic pairs with low relative kinetic en- 



2 



Chr. Samuelis el al: Optical manipulation of long-range interactions at the 3s+3p asymptote of Na2 



ergy. We extract quantitative information on light induced 
coupling from the induced level shifts. A similar laser con- 
trolled manipulation of the asymptotic ground state po- 
tential can be used to investigate schemes to modify the 
ground state s-wave scattering length a s . The advantage 
of the experiment at the 3s!/ 2 +3pi/2 asymptote presented 
below is the simpler detection via direct laser-induced flu- 
orescence, which is not possible in the case of the ground 
state asymptote. 

This article is organized as follows: First, an outline 
of our experimental setup and the excitation scheme is 
given (Section EJ. In the following Section |31 we focus 
on the extraction of quantitative data from the observed 
line profiles. Experimental results are presented in Sec- 
tionEI A theoretical model that simulates the investigated 
light induced energy shifts and line broadenings utiliz- 
ing a coupled-channel calculation including six channels 
is developed in Section in the framework of the Mapped 
Fourier Grid Hamiltonian (MFGH) method involving a 
complex potential |19U20j . In the concluding section we 
compare the results of the simulation with the experi- 
ment and give also prospects for further experiments and 
extended simulations. 



2 Experiment 

We briefly review the experimental setup and focus on 
the changes made in comparison to |21| were we studied 
asymptotic sodium A state levels within our molecular 
beam experiment. For the investigation of light induced 
energy shifts below the asymptote of the A state, we start 
from molecules in a well collimated beam (molecular ve- 
locity of « 1000 m/s) which are initially in the lowest vi- 
brational levels (vx — 0, 1) of the AT 1 X'+ ground state (see 
Figure GJ. In a first interaction zone (Figure EJ) a Franck- 
Condon pumping step is applied with laser LI to cre- 
ate population in thermally unpopulated vibrational levels 
vx ~ 31 of the ground state. For this purpose a dye laser 
(sulforhodamine b) operating at about 615 nm is used. Its 
frequency is stabilized on the maximum of fluorescence of 
the selected molecular transition. 

Starting from such vibrationally excited ground state 
levels the A state asymptote is reached with a tunable 
dye laser L2 at 532 nm operating with coumarin 6. It is 
applied in a second interaction zone about 0.35 m down- 
stream from the first zone. The fluorescence from the ex- 
cited levels is monitored by a photomultiplier. An OG 570 
color glass filter suppresses the scattered laser light from 
L2. 

For the manipulation of bound levels close to the A state 
asymptote a third laser L3 is used. It is an extended cav- 
ity diode laser in Littrow configuration operating near the 
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3 2 D 3 / 2 transition at 818 nm. 



It 



atomic sodium 3 2 Pi 

is overlapped with L2 by a dichroitic beam splitter and fo- 
cused into the second interaction zone. In that way a cou- 
pling of the A X E+ state to the 4 1 E+ and the 2 1 L T g state 
at the 3s+3d asymptote is created. The waist radius is 
ss 80/iin for L2 and about w 100/im for the laser L3. Over- 
lap and diameter of both beams are controlled with a CCD 
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Fig. 1. Potential scheme of Na2. Laser LI transfers popula- 
tion to a desired vibrational level (vx = 31), so that the 3s+3p 
asymptote can be reached via the second laser L2. The oblique 
arrow indicates that this laser is scanned during the experi- 
ment. A coupling of asymptotic A state levels to the 3s+3p 
asymptote is induced via the coupling laser L3. The fluores- 
cence from the A state is detected with a lock-in detector syn- 
chronized to the chopper frequency of laser L3. 
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Fig. 2. Experimental Setup. A well-collimated molecular beam 
is crossed by three laser beams. Laser LI is stabilized on the 
molecular fluorescence in the first interaction zone, while the 
fluorescence in the second interaction zone is recorded. The 
perturbing laser L3 is locked on a stabilized cavity. The signal 
from the second interaction zone is processed with a channel 
switch to record a manipulated (L3 on) and unmanipulated 
(L3 off) spectrum, simultaneously. 



camera. Both laser beams are linearly polarized, and their 
polarization axes are parallel to the molecular beam. 

For the determination of the atomic transition fre- 
quency 3 2 Piy 2 — ► 3 2 D 3 / 2 we apply optogalvanic spectros- 
copy in a conventional sodium hollow cathode lamp as 
reference. By this means a Doppler broadened spectrum 
of the atomic transition can be obtained. The position of 
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Fig. 3. A spectrum w 2.5 GHz below the asymptote of the 
A state. LI was stabilized on the R(6)(17-0) transition, and 
L2 is tuned across the P(6)/R(6) transitions vx — 31 — * va — 
178, 179. The dotted line is a signal with the coupling laser 
switched off, while the solid one corresponds to the situation 
when the coupling laser is on (Jl3 ~ 74 W/cm 2 ). The back- 
ground results from scattered light of L2. The additional off- 
set in the solid spectrum is due to additional scattered light 
from the coupling laser. This laser is stabilized in frequency 
with a detuning of 75 MHz to blue with respect to the atomic 
sodium 3 2 Pi/2 ~+ 3 2 D 3 / 2 transition. For both vibrational lev- 
els the line profiles of the P and R line are asymmetric. Each 
line has a shoulder on the low frequency (left) side originating 
from underlying hyperfine structure. 



the peak can be determined with an uncertainty of about 
30 MHz. For a good long-term frequency stability of L3 
it is stabilized to a 150 MHz marker cavity that itself is 
locked to an iodine stabilized HeNe laser [221 • L° n g term 
drifts of L3 less than 1 MHz/h are achieved in this way. 
The cavity is locked to the HeNe laser using an acousto- 
optic modulator (AOM) which allows interpolation be- 
tween the 150 MHz markers. Hence, we are able to deter- 
mine the relative detuning with an absolute uncertainty 
in the order of 1 MHz as long as the cavity is locked to 
the HeNe laser (see Section 0J. 

The effect of L3 on the energetic position of asymptotic 
levels in the A state is observed by a switching technique: 
L3 is modulated with a mechanical chopper at a frequency 
of about 1 kHz. The current of the photomultiplier in the 
second interaction zone is amplified with a fast current 
amplifier and then fed into an electronic switch, triggered 
by the chopper signal. It divides the signal into two out- 
put channels: The first channel is correlated to 'L3 off', 
the second to 'L3 on' if the switching phase is adjusted 
properly. 

A typical recorded spectrum, tuning laser L2 while 
keeping LI and L3 fixed, is shown in Figure OH Several 
scans are averaged to improve the signal-to-noise ratio of 
the recording. For the presented scan the total integration 
time per data point was Is. The Franck-Condon pumping 
laser LI is usually stabilized to the A-X R(6) (17-0) tran- 
sition. The population in va =17, J a =7 decays to Jx =6 
and 8 preferentially of vibrational levels around vx =31. 



Starting from ux=31, J x =6, L2 stimulates the transitions 
P(6) and R(6) to the asymptotic vibrational levels va — ^78 
and 179. The spectrum for every va thus consists always 
of two angular momentum components J a — Jx ± 1 for 
each va, to which we refer as P line [J a = Jx — 1) and 
R line (J a = Jx + 1) in the following. The spectrum is 
shifted to red (lower frequency) by laser L3, the lines are 
broader and show asymmetry. 



3 Extraction of Line Shifts and Line 
Broadenings 

In this section we explain how we extract information on 
the line shift A between the manipulated and the unma- 
nipulated trace in a spectrum like in Figure We deter- 
mine the induced light shift by applying a line shape fit. 
We deal with the two channels — with and without the 
coupling laser L3 — simultaneously because the two traces 
can be described by several common line shape parameters 
and few additional parameters describing the line shift. It 
is useful to concentrate on simultaneous simulation of P /R 
doublets because the two lines overlap for large induced 
line shifts. The line broadening is determined within the 
same line shape fit. 

For the description of the spectrum several hyperfine 
components underneath each rotational line have to be 
taken into account: Due to nuclear spin statistics both P 
and R line for an even Jx have contributions from nu- 
clear spins I~0,2. For the latter different orientations 
with respect to J a appear, which results to 1+5 hyper- 
fine components. For the electric dipole transitions I and 
Mi remain unchanged to good approximation. Each con- 
tributing hyperfine component is described by a standard 
profile 



A[A ,,,, ] = TT ^ 



A 



with 



X 



+ (1 — j3 — a)X 4 + f3 6 X 6 



v - vq 
HWHM' 



(1) 



(2) 



The parameters a, (3 and HWHM (half width at half 
maximum) are equal for all lines contributing to the un- 
shifted spectrum, a and (3 are line profile parameters, 
which depend on the contribution of the Lorentzian (a = 
1, (3 = 0) and residual Gaussian shapes to the line shape. 
They are determined by the collimation ratio of the ex- 
perimental apparatus (1:1000) and possibly by frequency 
jitter. For the profile fits presented here, a has been kept 
fixed at 1 while (3 is variable. 

The fit results show that the observed line profiles can 
be simulated using the scalar and the tensorial nuclear 
spin-spin interactions. The center frequency vq and the 
amplitude of each line component thus does not only de- 
pend on the rotational quantum number J a and the nu- 
clear spin /, but also on the total molecular spin F = J+I. 
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The total contribution diagonal in J of the hyperfme in- 
teraction to the line position is |23| : 



^hfs =(n = oji F\H BFS \n = ojif) 

I(I+l)-2i(i + l) 



4" 



+ d i(i + l)(2i + l)(-l) /+F+1 

r i 1/2 

x 30(27+1) (27 +1) 



(3) 



x (2J+ 1) 



J J 2 





1=1 =ii+i2 is the total nuclear spin, i = ii=i2 = 3/2 
is the nuclear spin of a single sodium atom and Wigner's 
symbols are used. The first term in equation © is the 
scalar part of the hyperfme interaction while the second 
one is the tensorial contribution of the hyperfme energy. 
The ^-parameter for the scalar and the d-parameter for the 
tensorial part depend on the states mixed by the hyper- 
fine interaction. At the 3s+3p asymptote a large variety 
of states is coupled to the A state, so that the functional 
dependence of the ratio of 6 and d does not have a simple 
analytic form. Thus, for the P/R line we let both parame- 
ters vary in the fit. The quadrupole and magnetic rotation 
hyperfme interactions are negligible in the present case of 
asymptotic levels. 

For the profile simulations, we assume that all hyper- 
fine components for a single J resulting from equation J3J 
are of equal intensity and thus we need two parameters 
for the description of the intensities of each P/R doublet. 
This assumption would be not sufficiently good for small 
angular momenta J, F. For the positions of P and R lines 
= v v + B v Ja(Ja + 1) + 7-hfs is used. For every va the 
unperturbed profile is composed by (for even Jx) 

U(v)=U oS + J2 E ( 4 ) 

J A = Jx±l 1=0,2 

Ja+I 

A[A (7,F, J A ),v,uo(I,F, J A )], 

F=\J A -I\ 

with U f[ to account for a background offset in the exper- 
imental trace. 

For the simulation of the recordings with the coupling 
light field on, we use the same parameters a and /3 but 
a different HWHM. It can be larger than without cou- 
pling laser due to predissociation broadening induced by 
the optical coupling to the upper states (4 1 i7 s h , 2 1 27 g , see 
discussion below). 

The offset of the trace recorded with coupling laser S s 
differs from the offset U g of the uncoupled trace as the 
detector sees some additional scattered light originating 
from the coupling laser. We introduce an additional pa- 
rameter for the description of the absolute intensities of 
both P and R line of the shifted spectrum. 

A coupling laser field causes the splitting of a line at 
vq into two Autler-Townes components at vq ± A. The 



line shift A and the relative weight of the Autler-Townes 
components of each line depend on the frequency of the 
coupling laser and on the Rabi frequencies of the coupling 
transitions. For P and R line the ratio of the intensities of 
the two Autler-Townes components are fitting parameters 
and are later called qj A . Due to reasons discussed below 
the shifted doublet is composed of various lines with dif- 
ferent shifts A. 

The geometries of the laser beams L2 an L3 need to be 
considered for a proper simulation: The detection proba- 
bility is proportional to the local intensity of L2, Il2{t), 
and as the dipole coupling is M dependent, the induced 
shift is a function of M and of the local intensity of L3, 
7 L3 (r). 

As the polarization vector of L2 and of L3 are defined 
in the laboratory frame and the transition moments of the 
molecules are given in the molecular frame a transforma- 
tion between those two systems has to be performed. We 
assume that the projection M of the angular momentum 
J onto to the laser polarization vector remains unchanged 
in the second interaction zone, and that all M are equally 
populated. The latter assumption is justified as we start 
from a thermal molecular beam and the light of LI is un- 
polarized. 

As P, R, and Q (state 2 1 77 g ) coupling contribute to the 
induced shift of a level for a specific M, the M-dependence 
is not the regular parabolic dependency as for a single P, 
Q, or R coupling like M 2 or J 2 — M 2 that results from 
the direction cosine matrices. We derived from a simula- 
tion with diabatic, rotationally corrected potentials and 
ab initio coupling strengths (see Section [5J that for large 
internuclear distances the shift A resulting from P, Q, and 
R coupling can be factorized into a function parabolic in 
M and linear in I^: 



Z}(M,7 3 )~7 L3 



1-6 



M 



(5) 



We derived b = 0.7 in our simulation. This result of poten- 
tial shift calculations at large internuclear distance is ap- 
plied to asymptotic levels because in such cases the main 
contribution of the overlap integral of the vibrational wave 
functions results from large internuclear separations. 

For the geometric model of the laser intensity distribu- 
tion we assume Gaussian TEMoo laser modes. Although 
the beams L2 and L3 are focused (waist w 2,3 ) into the sec- 
ond interaction zone, we can assume that we have plane 
waves because the Rayleigh ranges z 2 ^ 3 are larger than the 
radius of the detection zone p (w <C p <C z l R ). With good 
approximation we can assume that the molecular density 
is constant in the observed zone. Taking the cylindrical 
symmetry of the laser fields into account, the description 
of the space dependent interaction can be reduced to one 
coordinate: the radial distance r from the common axis of 
the laser beams L2 and L3. We define as an approximation 
to a Gaussian profile a distribution of a discrete number 
Pmax of equi-thick hollow cylinders. In each of them we 
assume a constant intensity 7£ 2 and 7^ 3 representing the 
average intensity of L2 and L3 in the cylinder p where p 
counts from the center to the outer range. 
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Fig. 4. Experimental and simulated spectrum of v A = 178, 
J a = 5, 7. The dotted trace is the unshifted experimental spec- 
trum. The dashed trace is the shifted experimental one (cou- 

- 3 2 D 3/2 resonance, 7 L 3 =74 



pling laser on the atomic 3 P1/2 



W/cm 2 ). The two solid traces correspond to the simulation of 
the spectra as explained in the text. 



With this model the profile of one Autler-Townes com- 
ponent of one rotational line Jx —> J A with an unper- 
turbed transition frequency vq can be described as follows: 



T* [v,u , J A, Aq] 



Jx 

E 

M=-Jx 



J A Ux 

-M M 



(6) 



E 

P =i 



2p-l 

Pmax 



A[v,v ±A(M,Il 3 ),A ] 



The Wigner 3j-symbol gives the relative detection proba- 
bility for a single M in the unsaturated case for L2 and 
'Ji~ is the statistical weight for the hollow cylinder p 

of molecules. Like for the unshifted spectra described by 
equation Q, the summation over all contributing lines 
needs to be performed for the shifted spectra, too: 

Ja + I 

s(u) = s oS + E E E ( ? ) 

J A = Jx±l 1=0,2 F=\J A -I\ 



q±T+ [u, v (I, F, J A ) ,J A ,A (7, F, J A ) } + 
(1 - g± ) T" [v, v Q (I, F, J A ) , J A , A (I, F, J A ) } 



The intensity ratios of the Autler-Townes components qj A 

are limited by < q^ A < 1 and are either almost or 
almost 1 because the detuning of L3 from the coupling 
transitions is rather large compared to the Rabi frequency 
(see above). 

In the following the derived parameter of the level shift 
corresponds to the maximal shift, i.e., the shift for M = 
molecules at the beam center, z4(0,7l3) in equation (JSJ. 

Figure 0] shows the result of a fit to an experimen- 
tal spectrum for v A — 178. The experimental traces are 
plotted in dotted (unshifted) and dashed style (shifted). 
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Fig. 5. Effect of the optical manipulation on the last vibra- 
tional levels of the A state. The coupling laser is 75(30) MHz 
blue detuned with respect the atomic 3 2 Pi/2 — >3 2 D3/ 2 tran- 
sition. The intensity is 43 W/cm 2 . The dotted spectrum is 
recorded in absence of the coupling laser while the solid one is 
recorded with coupling laser switched on. 



L3 was resonant to the 3 P 



1/2 



>3 2 D 



3/2 



transition and 
caused an intensity of Jl3 = 74 W/cm 2 . The criterion 
for the quality of the fit is the sum of squared difference 
between experimental and fitted points. The agreement 
between experimental and fitted line profiles is within the 
noise interval. The line shift obtained from the fit is -28.9 
MHz for J A =5 and -29.7 MHz for J A —7 in the above case. 
The uncertainty of the shifts can be estimated by varying 
the starting conditions of the fit and is determined to be 
less than 3 MHz here. 



4 Experimental Results 

The dependence of the light induced energy shifts of v A — 
173, . . . , 182, J A =5,7 on power and detuning of the cou- 
pling laser L3 has been investigated. First we will focus 
on the effect on the last bound levels directly below the 
A state asymptote. Figure [S] shows a spectrum of the lev- 
els v A — 181 to 184. The trace plotted in dotted style 
is recorded without the coupling laser while the second 
trace, plotted in solid style, is recorded in presence of the 
coupling laser field. The intensity of the coupling laser 
was set to 7l3 = 43 W/cm 2 and the frequency was stabi- 
lized 75(30) MHz blue detuned with respect to the atomic 
3 2 Pi/2 ^3 2 D 3 / 2 transition. All levels are shifted to lower 
frequencies or energies, note the position of the poten- 
tial asymptote. Additionally, the line profiles of the levels 
v A = 182 to 184 have remarkably changed: The lines are 
broader than without the coupling laser. The P and the R 
line start to overlap. The dissociation continuum has been 
shifted, v A = 184 disappears. 

In FigurcEJthe variation of the observed shifts of v A =178, 
J A —5 as a function of the intensity of the coupling laser is 
shown. The laser frequency of L3 was kept fixed 75(30) MHz 
blue detuned as before, but the intensity was varied from 
15 to 67 W/cm 2 . This corresponds to energy shifts in the 
central region of the beam up to 26.9(30) MHz. The de- 
pendence of the shift on the laser intensity is linear within 
the experimental accuracy. 
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Fig. 6. Dependence of the light shift of va =178, Ja =5 on 
the intensity of the coupling laser. Without coupling laser this 
level is bound by 2646 MHz. The frequency of the coupling 
laser is stabilized 75(30) MHz blue detuned from the atomic 
3 2 Pi/2 -^3 2 D 3 / 2 transition. The absolute uncertainty of the 
coupling laser intensity is 50% and originates from the resolu- 
tion of the beam view camera utilized for the determination 
the laser beam diameters. The relative uncertainty between 
different measurements is less than 5%. 



The error limits for the absolute intensity result mainly 
from the uncertainty in the determination of the focus di- 
ameter of the coupling laser beam. We use for this pur- 
pose a CCD camera with a pixel size of 23/im x 27/im, on 
which we project a 10% reflex of each laser beam of the 
second interaction zone. The CCD and the second interac- 
tion zone are equi-distant from the beam splitter, so that 
the beam diameter measured with the CCD is the same 
as the focus diameter in the second interaction zone and 
can be derived with an absolute uncertainty of one pixel 
size of the CCD chip which results in the 50% uncertainty 
for the determination of the coupling laser intensity. The 
relative change of the intensity in Figure is known with 
an uncertainty of 5%, limited by the uncertainty of the 
power meter. 

As an example for the dependence of level shifts on 
the detuning of the coupling laser experimental shifts of 
va —179, Ja —7 are plotted in Figure ITUl The reference 
frequency for the detuning is the energy difference 5u> from 
va~179, J a— 7 to the 3s!/2+3d3/2 asymptote. This is rea- 
sonable because in the region of the outer turning point of 
the A state vibrational wave function (where the coupling 
is induced) the potential of the upper 3s+3d asymptote 
has almost reached its asymptotic value. This is due to the 
1 / independence from a quadrupole-quadrupole coupling 
of a s-electron and a d-electron and is in contrast to the 
1 / 7? 3 -behavior of the resonant dipole-dipole interaction of 
the A state potential correlated to the s+p asymptote. For 
the series of measurements in Figure HU1 the laser intensity 
of the coupling laser was held fixed at 74 W/cm 2 . The 
circles in Figure HUI are results from simulations, which we 
will describe in the following section. 
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Fig. 7. Potential scheme of the dipole coupling of A 1 state 
to the 4 J I7+ and the 2 1 77 9 states. Dipole selection rules allow 
P and R coupling to the 4 1 £'j~-state, and P, Q, and R coupling 
to the 2 1 77 9 -state. The inset shows the rotational barriers at 
the 3s + 3d asymptote for J a = 7. 



5 Theoretical Interpretation 

The observed level shifts originate from the laser-induced 
dipole coupling of the A 1 27+ state to the 4 1 Z ,+ and the 
2 1 77 g states, both correlated to the 3s+3d asymptote (Fig- 
ure P). We use a dressed molecule picture similar to [21] 
for the simulation of the coupled system. After eliminating 
the explicit time dependence from the theoretical model, 
a mapped Fourier grid method |19j is used to determine 
line positions and light induced line broadenings. 

For the correct description of the dipole coupling we 
have to take care of the e/f symmetries of the involved 
rotational levels and the selection rules for an optically 
induced dipole coupling. Thus, the A state can only cou- 
ple to the 4 1 Z'+ state on a P and R line, while there is P, 
Q, and R coupling to the 2 1 7T g -state. For a proper sym- 
metrization, the latter state usually is expressed as a linear 
combination ^ITg) = -j= (\f2 = 1) ± (-1) J |^ = -1)). The 
plus sign gives the e component of the 77-state that can 
only be reached by P and R couplings, while the minus 
is the f component which couples to the A state by a Q 
line. For each M (projection of J onto space fixed axis 
neglecting hyperfine structure), a coupled system can be 
calculated. Thus, the system to be solved consists of six 
(except five for M=0, where the Q coupling for £ — > 77 
vanishes) coupled channels, as illustrated in Figure 

The coupled system is described by a six component 
wave function. For a simple notation we use G, S, and 
77 instead of A 1 ^, 4 1 Z ,+ , and 2 1 77 9 for the electronic 
part of the wave functions. With Born-Oppcnhcimcr basis 
functions, the coupled channel wave function thus can be 
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expressed as 



*(J2,t) 



&b J1 (R,r k ,t 



^~ 1] {R : v kl t 
^ ( iP(R,r k ,t 
\* { n J+1) (R,r k ,t 

f±XG(R,t) 4>g 
\xs(R,t) <j>E 
\xs(R,t) 0z 
RXn(R,t) <pn 
jiXn{R,t) <p n 

\jiXn{R,t) <p n 



J 

ThlR)\tj) \ 

r k ;R) lO-i) 
r k ;R) lO+i) 
r k ;R) 
r k ;R) IO) 

v k ;R) \Z J+1 )J 



where R is the distance of the two nuclei and rfc are the 
coordinates of the electrons in the molecular frame. 
are the basis functions for the rotation for an angular mo- 
mentum J. 0c(rj;i?), (j>s(ri;R), and (j)n{^i]R) are the 
electronic wave functions corresponding to the potential 
curves Ug, Us, and Un- The functions \G, Xs, and xn 
describe the radial motion in these potentials. 

The Hamiltonian of the coupled system divides into 
two parts: 



- mol 

H = H 



W. 



(9) 



r mol 



The molecular part of the Hamiltonian H describes 
the field free case, whereas the second part W describes 
the dipole coupling between the different electronic states. 
Both matrices are set up as a matrix of kxk blocks where k 
is the number of coupled channels. The diagonal blocks of 

H are the matrixes H™ 1 (1 < i < k), all non-diagonal 
blocks vanish. W is set up of Wy (1 < i,j< k). Wij 
represents the dipole coupling between the states of chan- 
nel i and j . W i:) - vanishes if i and j are equal or if none of 
the channels i,j is G. Operators denoted with underscore, 
i.e., H, represent operators that act on all channels (in our 
case five or six for all internuclear distances). In contrast 
operators without underscore represent operators that ei- 
ther act on one channel only or describe the coupling of 
exactly two channels (H™°' resp. Wy). 



H 



mol 



2fiR 2 



(10) 



is the molecular part of the Hamiltonian. Tj is the kinetic 
energy operator and Uj the electronic potential energy 
operator of the electronic state i. The third term is the 
rotational energy of the two nuclei relative to each other 
represented in Hund's coupling case (a). \i is the reduced 
mass of the molecule. 



The coupling laser L3 is assumed to be a plane and 
monochromatic wave with frequency ujl- Thus, in dipole 
approximation 



w (*)i 



D 



i,j(R) ■ £lE - (e 



-iuj L t\ 



(11) 



D 



t .j(R) is the dipole moment for the transition between 



(g) the electronic states <f>% and <pj resulting from the integra- 
tion over rfc according to equation ©. Di.j(R) vanishes 
for transitions between S- and 77-state because g <-> g 
couplings are forbidden by electric dipole selection rules. 
6l describes the unity vector in the direction of the laser 
polarization and Eq is the electric field strength. 

With this Hamiltonian, the time-dependent Schrodinger 
equation for the six-components wave function SS?(t) is: 



ih^(R,t) 



- mol * , , 

H +W{t) H?(R,t) 



(12) 




200 300 
internuclear distance R [A] 

Fig. 8. Asymptotic region of the coupled electronic potentials. 
The A state potential is dressed by a photon energy of the cou- 
pling laser L3. The energy difference Aui is due to the detuning 
of the coupling laser relative to the atomic transition. Slo is the 
relative position of one individual vibrational level va to the 
3si/2+3d3/2 asymptote within the dressed state picture. 



The radial wave functions for the nuclear motion in 
the electronic channels at the lower and upper asymptote 
are transformed via 



ti 2 [ J ( J + 1) + S (S + 1) + L (L + 1) - Q 2 - S 2 - A 2 ] . 



Xc{R,t) = XG(R,t)exp(iAojt), 
Xe{R, t) = xs(R, t) exp(-iw t), 
Xn(R, t) = xn(R, t) exp(-iuj t), 



(13) 



where ujq (see Figure d) is the transition frequency of the 
atomic 3 2 piy 2 ~* 3 2 d 3 / 2 -transition, Alo is the detuning of 
the coupling laser with respect to ujq (Ad = ojl — ujq). 
By the transformation l|13|) we shift (dress) the origin of 
energy for the electronic potentials, 



U G {R) = U G (R) + HAuj, 
U S (R) - U S {R) - huo, 
U n (R) = U n (R)-hLu . 



(14) 
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Hence, all three electronic potentials have almost the same 
asymptotic energies. In Figure |H1 the resulting molecular 
potentials dressed by one photon of a blue detuned (rela- 
tive to the atomic transition) coupling laser are plotted. 

Adapting the dressed potentials in the H™ , the Schro- 
dinger equation l|12|) for the vibrational motion in the six 
coupled channels can be written as 



continua. Here we recall only the main steps of the MFGH 
method: 

— In the standard Fourier Grid Hamiltonian (FGH) method, 
the total Hamiltonian for the coupled channels is ex- 



pressed in a basis on N plane waves, exp 



1T1 



at 



X J n \R,t) 
X J n(R,t) 
\x J n + \R,t)J 



( Kg* 1 

nn° n 
\hnjj 1 



nn^ 1 
o 







H 



mol 







o 
o 

xjmol 

IIj-i 



o 
o 
o 

umol 



o 







H 



-J 



x J G (R,ty 

x J z\R,t) 
xi +1 (R,t) 

xirHR,t) 

X J n(R,t) 
\x J n +1 (R,t). 

(15) 



where J is short for J a used in the other sections. 

To solve this multi-channel problem the explicit time 
dependence of the Hamilton operator can be eliminated 
if the rotating wave approximation |25| is applied. It is 
based on the suppression of the fast oscillating terms in the 
coupling matrix elements, which average out over times 
given by the vibrational period of the molecule. 

The coupling matrix elements between the G (origi- 
nally the A state) and electronically higher excited states 
S and II are expressed by time independent Rabi frequen- 
cies nP(R)/2n (j = U,n) where j3 = 0,±1 indicates the 
difference of the rotational quantum number between the 
coupled states. They are related to the electronic transi- 
tion dipole moment by 



1 E \D G ,j(R)\ ( J+/n 1 



2 

I / 2J L3 
ceo 



~M M 



\VgAR)\ [ J -mom 



(16) 



k = —(y — 1), . . . , 0, . . . , y, where L is the range of 
internuclear distances under investigation. The diago- 
nalization yields eigenfunctions represented as an ex- 
pansion over the same basis, the coefficients being the 
value of the wave function at each point of the grid of 
length L with TV equi-distant points. 
As the local de Broglie wavelength of wave functions 
varies by orders of magnitude between the small and 
\ the large internuclear distance range, an adaptative co- 
ordinate can be defined in the FGH framework, which 
maps the variation of the local kinetic energy. The 
j number of required grid points is reduced substan- 
tially, allowing accurate calculations of energy levels 
with large elongations, i.e., close to the dissociation 
limit. 

— The width of the energy levels interacting with a disso- 
ciation continuum is obtained after the diagonalization 
of the total Hamiltonian, including a pure imaginary 
potential (or an optical potential) V op t placed at large 
distances, which ensures absorbing boundary condi- 
tions for the outgoing waves. In a stationary approach, 
the dissociating wave functions behave asymptotically 
as Siegert states [20] ■ We choose the form proposed in 



E2: 



V'opt = A opt N opt exp 



-2L 



opt 



R — it. 



opt 



(17) 



where N ovi is a normalization factor, i? Q pt, 



I-opt and 



Here we take into account the earlier discussed ^/-depen- 
dence by the linear polarized light in the experiment. For 
each projection of M on the laboratory frame the eigen- 
value problem of the coupled channel system gives the sta- 
tionary solution and can be solved for each M separately. 
Levels with different M do not couple. 

The light induced energy shifts (the difference of cor- 
responding eigenvalues with and without light field) and 
line broadenings are computed using the Mapped Fourier 
Grid Hamiltonian (MFGH) representation, described in 
detail in The MFGH method has been designed in 
the context of photoassociation of cold atoms and cold 
molecule formation, where interactions at large internu- 
clear distances play the key role. In subsequent papers 
|2"oTl2Tj] . a complex potential has been introduced to treat 
the interaction of bound levels embedded into dissociation 



A opt are respectively the position, the length, and the 
amplitude of the optical potential. The values of these 
parameters are chosen according to the recommenda- 
tion of |27j . The imaginary part of the eigenvalues re- 
sulting from the diagonalization of the Hamiltonian 
matrix provides the predissociation line width 

In our calculations, we used a grid with L=1100 a.u.(l a.u. 
=0.5291772083 x 10~ 10 m), with typically 700 grid points, 
which induces a diagonalization of 4200x4200 squared 
matrix, using standard LAPACK routines. As we are look- 
ing for very small energy shifts, a high accuracy is needed 
for the eigenvalues, which is checked by studying their 
convergence depending on the contraction (equal to 0.2 
here) of the grid step. A contraction of 1 corresponds to a 
local grid density of four points per local de Broglie wave- 
length. The optical potential is set up to absorb outgoing 
waves with kinetic energies corresponding to the detuning 
Aio of the coupling laser L3, following the rules described 
in |27]. The parameter values are A opt — 5.6 x 10~ 7 a.u., 
L opt = 210 a.u., R opt = L - L opt , and A^pt = 13.22. 
Changes of ^4 pt, L and L opt in reasonable ranges do not 
affect the simulation results for the line positions and the 
line widths, so that their choice is appropriate to represent 
the loss of molecules by predissociation. 
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■ shift J =7 M=0 

• light induced prcdissociation width M=Q 
I L =75 W/cm 2 
Aaj=0.0375 an' 1 



Fig. 9. Simulation of light shifts and line broadenings of the 
last 25 vibrational levels of the A X S^ state calculated for Ja~7 
and an intensity of the coupling laser of 75 W/cm 2 . The cou- 
pling laser is assumed to be 0.0375 cm -1 blue detuned to the 
atomic 3 2 pi/2 — > 3 2 d3/2 transition. All levels of the A state are 
shifted to lower energies. 



For the simulations, we applied the following poten- 
tials: For the A 1 state a RKR potential (Rydberg- 
Klein-Rees) generated from experimental values is used. 
It is based on energies from The potentials for 4 1 Z'+ 
and 2 1 n g state for internuclear distances up to 35 a.u. 
are taken from ab initio calculations in |28| (method b) 
therein). The long range part is calculated from the dis- 
persion coefficients given in and is connected contin- 
uously differentiable to the inner part. 

The dipole moments used for the light induced cou- 
pling originate from calculations in |28| . They converge 
for internuclear distances above 100 a.u. to the atomic 
values. For the inner part below 14 a.u., we generated ad- 
ditional values of the dipole moment by smooth continu- 
ation of the ab initio curves. This inner part of the dipole 
moments has only minor influence on the results of the 
numerical simulation as the contribution to the overlap 
integral of the A state and upper states wave functions at 
small internuclear distances is almost negligible. 

For the simulation of line shifts of asymptotic A state 
levels, we compare the theoretical energy positions of the 
single channel calculation of the A state with the cor- 
responding values of the six channel system. A typical 
result of a simulation is presented in Figure [5J We dis- 
play the energy shifts (squares) of vibrational levels va — 
165, . . . , 184 with J a— 7 for a peak intensity of the coupling 
laser 7l3 = 75 W/cm 2 and a blue detuning of 0.0375 cm -1 

3 2 D 3 / 2 transition. 



with respect to the atomic 3 P1/2 

Moreover, the figure contains light induced line broaden- 
ings (circles) obtained from the simulations with imagi- 
nary absorbing potential. With the blue detuning of 0.0375 
cm^ 1 all levels bound by less than this energy starting 
with va — 181 show considerable predissociation broad- 
ening of up to 49 MHz. 

Similar calculations have been performed for all stud- 
ied laser intensitites Jl3 and detunings Aui of the coupling- 
laser. In Figure Hill the circles show a simulation of the de- 



-40- 

N 
X 

s 

= -60 
-80 



-100 



v=179J =7 

A A 

■ experiment ' L3 =74 W/cm 2 
O simulation /, =75 W/cm 2 



-4000 3000 -2000 -1000 

detuning of L3 8co [MHz] 



Fig. 10. Dependence of the light shift of va =179, Ja =7 on 
the detuning of the coupling laser relative to a hypothetical 
transition (va =179, Ja~7) — > 3s 1 / 2 +3d 3 / 2 asymptote (see also 
Figure |7| and Figure |HJ . The squares are experimental points 
taken for a peak intensity of 7l3 =74 W/cm 2 . The circles are 
from a simulation with 7l3 =75 W/cm 2 . 



pendence of level shifts for va =179, Ja —7 on the detun- 
ing Slo of the coupling laser with respect to a hypotheti- 
cal transition (va =179, Ja =7)^3s 1/ / 2 +3d 3 /2 asymptote. 
Experimental points are plotted by squares. 

The agreement between experimental and simulated 
line shifts is convincing, but one has to keep in mind that 
due to the uncertainty in the determination of the cou- 
pling laser beam waist, we have a 50% uncertainty in the 
laser intensities applied for the simulation. For detuning 
less than 1000 MHz our fitting procedure fails to deter- 
mine line shifts from the experimental spectra. The exper- 
imental signal of the shifted trace decreases and broadens 
remarkably towards the resonance frequency. A typical 
spectrum showing this effect for va = 179, J a = 5, 7 is 
presented in Figure ITT1 The theoretical model applied for 
the fits of line profiles is no longer able to simulate the 
observed line profiles properly. 

The deviations in Figure ^2 are not noise; the slow 
variation is reproducible. For a proper simulation, the fit 
should include the fine and hyperfine structure on the op- 
tical coupling. 

A comparison of the line broadening observed in the 
experiment with theoretically predicted line broadening 
due to laser induced predissociation is shown in FigurelT2l 
The experimental broadenings are calculated from the fit- 
ted FWHM in presence of the coupling laser minus the 
FWHM without light induced broadening. In the experi- 
ment we observe a significant line broadening, which has 
a maximum if the coupling laser is on resonance with a 
hypothetical transition frequency from the level under in- 
vestigation to the 3s!/ 2 +3d3/2 asymptote (Su — 0). The 
broadening effect for Scu < cannot be described by the 
laser induced predissociation calculated with our six chan- 
nel model because the dissociation channel is not yet open. 
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Fig. 11. Example of an experimental scan of va =179, J a =5,7 
when the coupling laser is on resonance (Slu = 0) with respect 
to an hypothetical transition va =179, J a =7— *3s 1 / 2 +3d 3 /2 
asymptote. The intensity of L3 was Il3=74 W/cm 2 . The dot- 
ted trace is taken without the coupling laser, while the dashed 
one is the manipulated spectrum. The solid traces represent 
fits on the experimental data with the model presented in Sec- 
tion [3] 
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Fig. 12. Dependence of the line broadening of va = 179, 
J a = 7 on the detuning Slu of the coupling laser (see Figure|HJl. 
The squares represent the difference of the two FWHM fitting 
parameters. The circles are from a simulation for the experi- 
mental peak intensity of 7l3 =74 W/cm 2 , neglecting the fine 
and hyperfine structure of the asymptote 3s+3d. The threshold 
is only given for the discussion in the text. 
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Fig. 13. Hyperfine potentials at the asymptote 3s + 3d for a 
total angular momentum F = 0. 



adiabatic asymptotic potentials are shown in Figure 1131 
From the A state, which correlates to the 3si/ 2 {f = l)+3pi/ 2 
asymptote, a dipole transition is only possible to the 
3si/2(/ = l)+3d 3 / 2 asymptote. Thus, the light induced 
shifts of an asymptotic A state level shows a frequency 
dependence which can be described by an effective cou- 
pling to only this asymptote. In the region above the 
3si/2(/ = l)+3d5/2 asymptpote, which is the lowest asymp- 
tote due to the inverted fine structure, a predissociation 
can appear. This will lead to an increased line width of the 
A state levels under investigation starting at a detuning 
Suj > — \Eps\- This threshold is included in Figure IT21 

A second effect leading to a predissociation and thus 
a decreased lifetime of the 4 1 Z'+ state is the vibrational 
coupling of this state to lower lying molecular states with 
the same symmetry. In Figure ITU the region of an avoided 
crossing with the 3 1 Z'+-state is marked by a dashed box. 
The state correlates to the 3s+4s asmyptote so that 

in the region above this asymptote a predissociation to 
this asymptote is possible. Therefore, the optical coupling 
of the A state to the 4 1 i7+-state can also induce a predis- 
sociation of A state levels to the 3s+4s asymptote. 

Due to limitations of the memory available on our com- 
puters the two additional predissociation channels can not 
be implemented in our coupled channel code. Hence, it is 
difficult to estimate the contribution of the two effects to 
the experimentally observed line width. 



6 Conclusion 



There are two different effects, which can lead to the 
experimentally observed line broadening for detuning Slu < 
0. The molecular structure of the states at the 3s+3d 
asymptote is not as simple as we have assumed for the 
simulation. Due to the fine structure energy of the 3d 
atom and hyperfine energy of the 3s atom, which have 
almost the same magnitude (_Eps(3d)=-0.0498 cm -1 
£hfs(3s)=-0.0590 cm" 1 EH]) the 3s+3d asymptote splits 
into four different hyperfine asymptotes. The calculated 



A high resolution molecular beam experiment has been 
used to investigate the effect of a light induced coupling 
on the last vibrational levels of the sodium A state. The 
coupling is induced by a laser field which is near res- 
onant to the atomic sodium 3 2 P]y2 - * 3 2 D 3 / 2 transi- 
tion and thus couples the A 1 !]^ state to the 4 1 I7+ and 
2 1 n g state, which are both correlated to the 3si/2+3d 3 / 2 
asymptote. From the experimental data we see that with 
the usual laser intensities available from lab size laser 
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Fig. 14. Molecular potentials with Xg~-symmetry. The shape 
of the 4 1 X , ^"-potential is due to an avoided crossing with the 
3 X^-state that correlates to the 3s + 4s asympote |28) . 

sources (~ 50 W/cm 2 ) light induced level shifts in a 
molecular system in the order of several 10 MHz can be 
induced. Furthermore, we observe an increased line width 
of the A state levels under investigation due to light in- 
duced predissociation. The number of bound states of the 
A state could be changed by the influence of the coupling 
laser by one or more units. Thus, in the picture of collid- 
ing atoms, the scattering phase was altered by more than 

7T. 

We developed a theoretical model of the light induced 
coupling, which describes the observed level shifts in good 
agreement with the experimental data. It is based on a 
six channel calculation using the Mapped Fourier Grid 
Method for the solution of the coupled channel eigenvalue 
problem. We use ab initio data for the electronic potentials 
at the 3s 1 / 2 +3d 3/ <2 asymptote and for the dipole moments 
for the description of the light induced coupling. For the 
simulation of the induced line shifts of deeply bound states 
it is not necessary to take into account the fine and hy- 
perfine structure of the 3s+3d asymptote. 

The interpretation of the light induced line broaden- 
ing effects needs a more detailed analysis of the molecular 
structure at the 3s+3d asymptote. The predissociation of 
the 4 1 Z'+ and 2 1 n g state due to fine structure, hyper- 
fine structure, and vibrational motion must be taken into 
account. 

This research is in progress with a similar experiment 
that couples the sodium ground state asymptote with the 



3s+3p asymptote [32]. Here we want to investigate sys- 
tematically a laser induced coupling between the X 1 i7+ 
ground state of sodium to levels in the A state. This ex- 
perimental situation is of very high interest, as in most 
cooling and trapping experiments lasers (either for a MOT 
or dipole trap) are present, which are slightly detuned 
from the atomic resonance. These experiments show very 
complex dynamical Stark effect which asks for further ex- 
tensions of our model system. 
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